zeta_values = [0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1];  % 创建包含不同阻尼比的向量zeta_values

t = 0:0.01:30;  % 创建一个从0到30，步长为0.01的时间向量t

step_responses = zeros(length(t), length(zeta_values));  % 创建一个与t和zeta_values大小相匹配的全零矩阵，用于存储步响应数据

for i = 1:length(zeta_values)  % 对于每个阻尼比值
    zeta = zeta_values(i);  % 从zeta_values中获取当前阻尼比值
    num = 1;  % 系统的分子多项式为1
    den = [1, 2*zeta, 1];  % 系统的分母多项式
    sys = tf(num, den);  % 创建传递函数模型
    step_responses(:, i) = step(sys, t);  % 计算系统在t时刻的阶跃响应，并存储到step_responses中的相应列中
end

figure  % 创建一个新的图形窗口
ribbon(t, step_responses)  % 使用ribbon函数绘制彩色带状图，表示不同阻尼比下的阶跃响应
title('归一化二阶系统阶跃响应')  % 设置图表标题
zlabel('y(\zeta,t)')  % 设置z轴标签
legend('\zeta = 0.1', '\zeta = 0.2', '\zeta = 0.3', '\zeta = 0.4', '\zeta = 0.5', '\zeta = 0.6', '\zeta = 0.8', '\zeta = 1')  % 设置图例
view([130,30])  % 设置视角
shading interp  % 设置渲染方式为插值
colormap("jet")  % 设置颜色映射为"jet"
light  % 添加光照
lighting phong  % 设置光照模型
box on  % 显示坐标轴框线